Absence of magnetic ordering in NiGa 2 S4 
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Triangular-layered NiGa2S4, contrary to intuitive expectation, does not form a noncollinear an- 
tiferromagnetic structure, as do isoelectronic NaCrC>2 and LiCr02. Instead, the local magnetic 
moments remain disordered down to the lowest measured temperature. To get more insight into 
this phenomenon, we have performed first principles calculations of the first, second end third neigh- 
bors exchange interactions, and found that the second neighbor exchange is negligible, while the 
first and the third neighbor exchanges are comparable and antiferromagnetic. Both are rapidly 
suppressed by the on-site Hubbard repulsion. 



NiGa2S4 occurs in the layered structure where the 
main motif is a triangular layer of Ni 2+ ions surrounded 
by edge-sharing S octahedra, forming a trilayer S-Ni-S 
with the rhombohedral stacking ABC, and the trilayers 
are separated by the gallium oxide layers. Ni 2+ has the 
electronic configuration t^e^, an d one expects it to be 
insulating and magnetic with S = 1, and the magnetic 
moment per Ni being somewhat less than 2 [is- Indeed, 
this is exactly what happens with structurally similar 
transition metal oxides with a transition metal in the 
d & configuration, such as NaCrC>2 or LiCr02. 3d metal 
ions in this confuration do not have an orbital moment, 
therefore one expects a vanishing single-site anisotropy 
and magnetic interactions reasonably well described by 
the Hciscnberg model. In the nearest neighbor approx- 
imation this leads to noncollinear ground states, with 
neighboring spins pointing roughly at 120° to each other. 
Indeed, this is what has been observed in the above- 
mentioned chromates. 

NiGa2S4, on the other hand, has attracted substan- 
tial recent interest exactly because the experiments in- 
dicate absense of any long-range magnetic ordering 1 . 
Several explanations have been proposed, such as full 
cancellation of the nearest neighbor exchange and frus- 
trated competition between the 2nd and the 3rd neigh- 
bor interactions^, or biquadratic exchange-. These, how- 
ever, impose very severe quantitative restriction on the 
exchange parameters, which seem quite unrealistic. 

In order to elucidate magnetic interactions in this sys- 
tem we have performed ab initio density functional the- 
ory (DFT) calculations of the electronic structure and 
magnetic energies of NiGa2S4 and found that indeed con- 
dition required by either explanation are very unlikely to 
be satisfied, however, the magnetic interactions are very 
rich and (as conjectured in Refi) long range, so that tak- 
ing into account three neighbor shells is indispensable. 

For the calculations, the experimental crystal 
structure^ was used. A full potential linear augmented 
plane wave code* was used with a gradient approxima- 
tion for exchange and correlation 5 .. The calculations, as 
expected, render an insulating band structure shown in 
Fig. 1. As one can see, Ni d(e g ) bands are fully po- 
larized, a small gap opens (as usual, the absolute value 
of the gap in the DFT cannot be taken very seriously), 
the magnetic moment per Ni is 2 \xb-, and in the calcu- 




FIG. 1: LDA band structure of the ferromagnetic NiGa2C>4. 
The Ni character is emphasized by the size of the circles. 



lations this moment, not unexpectedly, resides entirely 
in the NiS2 layer. Interestingly, more than 20% of the 
total magnetic moment is located on the S sites. This 
creates a ferromagnetic interaction between the nearest 
neighbor nickels, which in the theory of strongly corre- 
lated magnetic systems is known as "ferromagneting 90° 
exchange" . Note that in the DFT the Hund rule energy 
is approximated as f dr/(r)m 2 (r)/4 rj N^IjM?/4, where 
m(r) is the total spin density, I and Mi are the Stoner 
factor and the total magnetization of the atom i, and this 
energy therefore does not depend on the Ni-S-Ni bond an- 
gle. Let us estimate this interaction^. Two Ni neighbors, 
when their spins are parallel, induce a magnetic moment 
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of ~ 0.2 \ib on each of the two bridging sulfurs, gaining 
an additional magnetic energy of 2/g0.2 2 /4 « 20 meV 
(The Stoner factor of the sulfur ion can be estimated as 
described in Refill and is about 1 eV). Defining J as half 
of the energy for flipping a bond (E nn , = J nn >S n • S„<), 
we find a ferromagnetic contribution to J of the order of 
10 meV. This is not a small energy, but it has to compete 
with also large conventional antiferromagnetic superex- 
change. The latter in DFT is 2t 2 /I Mi, where t is the 
effective Ni-Ni d — d hopping and Ijyi is the Hund en- 
ergy cost of exciting an electron with the opposite spin, 
In ~ 0.8 eV. However, taking into account on-site Mott- 
Habbard correlation mandates substiting the Stoner I 
in this expression by the Hubbard U, which is at least 
4 times larger. While Ni in NiGa2S4 is not necessarily 
strongly correlated, there is no doubt that the energy of 
adding an electron is substantially underestimated in the 
DFT, and hence the AFM superexchange overestimated. 
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FIG. 3: (color online) Calculated exchange constants for the 
first three neighbor shells in NiGa2 04, in meV, as a function 
of the Hubbard U, assuming the intraatomic J =0.07 Ry. The 
additional entry at U = corresponds to J = 0. 




FIG. 2: (color online) Different magnetic patterns used in 
calculating the exchange constants. Squares (circles) indicate 
up (down) moments within in the supercell, thin dark (thick 
light) line antiferro- (ferro-)magnetic bonds. The first pattern 
corresponds to the v3 x -^3 supercell indicated by the dash 
lines, the second, third and fourth patterns to a 2 x 1, 3 x 1, 
and 2x2 supercells, respectively 

Let us now investigate numerically the magnetic inter- 
actions in NiGa2S4. First, we want to make sure that our 
conjecture about the absence of magnetic anisotropy is 
indeed correct. This can be adressed by running fully 
relativistic calculations imposing different magnetic field 
directions and comparing energies. The result is that 
both energies differ by at most 0.03 meV. That is to 
say, the single site anisotropy is not an important fac- 
tor in NiGa2S4. Having established that, we have com- 
puted several different collinear magnetic structures, as 
shown in Fig. 2. R mapped onto a Heisenberg model 
with three nearest neighbor interactions, these give the 
exchange constants of 8.4, 0.3 and 4.1 meV (defined so 
that the total energy is equal to sum over all bands of 
JijSiSj^ \Si\ — 1) for the first, second and third neigh- 
bors, respectively. 

Several observations are in place. First, in LDA, while 



the second neighbor exchange is negligible, the third one 
is sizeable and comparable with the nearest neighbor 
exchange (it was suggested already in Ref. 0] that the 
3rd neighbor exchange may be anomalously large in this 
compound). This can be traced down to anomalously 
large 3rd neighbor hopping, which is in fact generic for 
triangular layers of transition metal oxides. Indeed, if 
the metal-oxygen bonds form precisely 90° angles, the 
strongest nearest neighbor hopping channel, e g — p — e g 
(or, in compounds like Na^CoC^, t2 g — p — t2 g ) is for- 
bidden by symmetry, however, a third neighbor path, 
e g — p — t2 g — p — e g is fully allowed and in fact has the 
most favorable geometry (FigdJ)^. This creates a possib- 
lity for a sizeable superexchange of the order of, as usu- 
ally, t^jf/A, where i e // is the effective hopping that ap- 
pears after all intermediate states are integrated out and 
A is the energy required to flip the spin of a metal ion. In 
LDA, t eff is of the order of t 3 pd J(E d - E p )/(E Eg - E t2g ), 
and A is of the order of the Stoner (Hund) parameter, 
< 1 eV. In the Hubbard model, on the other hand, A is 
set by the scale of Hubbard U ~ 4 — 6 eV. As usually, 
the real life is somewhat in between, meaning that the 
exchange constants in LDA are likely overestimated. 

This can be easily demonstrated using the LDA+U 
method that takes into account the Mott-Hubbard cor- 
relations in some crude approximation (Fig. [3]). To get 
an idea of the overall scale of the picture we have esti- 
mated the value of U using the quasiatomic loop in a 
standard Linear Muffin Tin Orbital package, as decribed 
in Ref. 7, and obtained U ~ 0.3Ry. This simplified 
method is known to underestimate U, therefore we have 
carried out calculations up to U = 0.45 Ry. At that max- 
imal value of U all three exchange constants practically 
vanish, within the accuracy of the calculation. In fact, 
the nearest neighbor constant at U = 0.45 Ry becomes 
negative (—0.04 meV), but for all practical purpose it 
may be considered zero. Interestingly, at this value of 
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FIG. 4: (color online) Different possible superexchange paths 
within MS2 trilayers. Thin dashed triangles indicate the tri- 
angular Ni and (in the left panel) S layers. Thick (red online) 
lines show the exchange paths. The solid parts correspond to 
short bonds (Ni-S, 2.422 Aand the shorter S-S bond, 3.212 
A), the dashed ones to the longer S-S bond, 3.626 
AA. 

U the sum of our calculated exchange constants over all 
(six) bonds gives 1.1*6=6.6 meV=76 K, to be compared 
with the Curie- Weiss temperature of 78±1 K— (the cal- 
culated sign is antiferromagnetic, in agreement with the 
experiment). 

What prevents the system from ordering remains un- 
clear. Nakatsuji et al^ conjectured that the ratio of 
the nearest and the third neighbor exchange (no second 
neighbors) is « —0.2, while our ratio, at large U, is es- 
sentially zero. To this point it should be mentioned that 
neither the accuracy of LDA+U functional is sufficient 
to make firm statements with a precision of a fraction 
of a meV, not the three-shell isotropic Heisenberg model 
is accurate to that extent. If the calculated numbers at 
U = 0.45 Ry are off by « 0.2 meV this would be enough 
to bring the calculated numbers in consistency with the 
Nakatsuji et aUs~ model. 

Finally, one should keep in mind that the actual or- 
dering temperature, if any, is suppressed by the 2D char- 
acter of magnetic interactions. Indeed, NiGa2C>4 has an 
extra Ga02 layer compared to typical layer oxides of the 
formula ABO2 and therefore the surexchange interaction 
between the layers is reduced. We have estimated this 
interaction by comparing the energies of the fully fer- 
romagnetic ordered structrure and the A-type antiferro- 
magnetic one (FM layers stacked antiferromagnetically) . 
The former is higher by about 1 meV (J =0.5 meV) in 
LDA and by about 0.3 meV (J = 0.15 meV) in LDA+U 
(U = 0.3 Ry). Importantly, this superexchange is addi- 
tive with respect to all possible hopping paths between 
the layers that include not only hopping from a Ni to the 
other Ni right above, but also to a large number of neigh- 
boring Ni sites in the next planed The real exchange 
coupling between the two antiferromagnetic planes will 
be additionally reduced. An estimate of J± ~ 0.05 meV 
seems reasonable. Of course, although a computer code 
calculates the numbers with arbitrary precision, the ac- 
tual physical approximations used in the calculations pre- 



clude statements about energy differences of the order of 
0.05 meV. For all practical purposes, it may be less than 
0.01 meV, which, of course would make long-range order- 
ing at the experimentaly probed temperatures (a fraction 
of a Kelvin) impossible. 

To conclude, we have calculated the three near- 
est neighbor exchange constants in NiGa204 by map- 
ping LDA+U calculations onto the isotropic Heisenberg 
model. We found that (a) an anomalously large third 
neighbor coupling exists in the system that can be traced 
down to superexchange via occupied t2 g orbitals, (b) 
upon including correlation effects in terms of Hubbard 
U, all exchange constants decrease rapidly, according to 
general superexchange intuition, but the first neighbor 
exchange is more rapidly suppressed than the third neigh- 
bor one, (c) at U = 0.45 Ry (6 eV) the nearest neighbor 
exchange is entirely suppressed (and we could not ex- 
clude, based on our calculations, that it does not become 
slightly ferromagnetic), while the third one is exactly the 
right magnitude to explain the observed Neel temper- 
ature. On the other hand, U = 6 eV, at least on an 
intuitive level, seems to be too large for Ni in such an 
enviroment. At more realistic f/'s, such as 3-4 eV, both 
interactions remains firmly antiferromagnetic. It should 
be kept in mind that energy differences on the order of 1 
meV are on the borderline of many approximations used 
in our analysis. It is possible that several weak effects 
conspire to prevent the system from ordering. First of 
all, in particular the nearest neighbor interaction results 
from cancellation of two considerably stronger superex- 
change interactions of the opposite signs: the AFM su- 
perexchange and the FM superexchange due to the Hund 
rule coupling on S. Albeit we see no obvious reason for 
the DFT to underestimate the latter, such a possibil- 
ity cannot be excluded. Second, we did not make any 
attempt to evaluate further exchange interactions be- 
yond the third shell. While there is no special mech- 
anism making them sizeable (as opposed to the third 
neighbor interaction), again we cannot prove that nu- 
merically. Third, while the calculated on-site anisotropy 
is very small, other effects beyond the isotropic Heisen- 
berg model, such as biquadratic exchange, dipole-dipole, 
multispin interactions, eto^ while small, may be not neg- 
ligible on the backgound of the strong cancellation of the 
FM and the AFM superexchange, and may possibly cre- 
ate additional frustration in the system. Finally, fourth, 
interplanar coupling between the antiferromagnetic non- 
collinearly correlated planes is at most a fraction of a 
Kelvin, and possibly even smaller. This additionally su- 
presses long range ordering. 

The author is grateful to D.I. Khomskii for useful dis- 
cussion and to C. Broholm for bringing up this problem. 
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